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Abstract. In the first paper of this series we explored the case where a quark-nova ejecta forms a degenerate shell, supported 
by the star's magnetic field. Herein, we consider the case where the ejecta has sufficient angular momentum to form a torus, 
and we show that the density and temperature of the torus are such that it will remain degenerate throughout it's lifetime. We 
go on to discuss the evolution of such a torus and apply it to AXPs, specifically 1E2259+586 and 4U0142-I-615. As it turns 
out, using our model we can account for many of the observations of these objects including the quiescent phase luminosity, 
and blackbody temperatures during both quiescence and bursting. Furthermore, for 1E2259+586 our model explains the steep 
and slow decay components seen in the burst lightcurve, as well as the rotation period glitches and enhanced spin-down rate. 
We also estimate the mass of the degenerate torus to be of the order of 10"* Mq, and speculate that the observed optical/infrared 
emission from 4U0142+615 might be a signature of the thin degenerate torus we describe here. 
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1. Introduction 

Discussed herein is the novel idea that a torus composed of de- 
generate matter from a quark-nova event could be responsible 
for features of Anomalous X-ray Pulsars (AXPs). As described 
in the first paper of this series (Ouyed, Leahy & Niebergal 
2006; hereafter referred to as OLNI), during a quark-nova the 
degenerate crust of a neutron star is blown off, leaving behind 
a quark star (QS) surrounded by left over, highly-metallic de- 
generate matter. 

In OLNI we discussed one possible fate of this matter; the 
case where the ejected crust had insufficient angular momen- 
tum to escape the QS's gravitational pull. Thus, it would either 
balance with the QS's magnetic field and form a co-rotating 
shell, or fall back entirely onto the QS. Using the shell sce- 
nario in OLNI, we were able to explain many features of AXPs, 
Soft Gamma-ray Repeaters (SGRs) and X-ray Dim Isolated 
Neutron stars (XDINs). Also, by appealing to the idea pre- 
sented in Niebergal et al. (2006), where a QS in the ground 
Color-Flavor Locked (CFLfl phase behaves as a type II su- 
perconductor and emits X-rays through magnetic vortex expul- 
sion, a good fit to X-ray luminosities of AXPs and SGRs was 
made. 

In this second paper the fate of the ejected neutron star crust 
is considered when the quark-nova compact remnant has ini- 
tial conditions (i.e. rotation period, magnetic field, and shell 
mass) such that the propeller mechanism acts. This, we argue. 



Send offprint requests to: ouyed@phas.ucalgary.ca 
' See OLNI on CFL and the assumptions in our model. 



results in a thin degenerate torus forming from the ejected mat- 
ter rather than a shell as in OLNI. As it is uncertain what initial 
conditions the quark-nova compact remnant will have, we feel 
that this second paper represents a more complete investigation 
of the fate of the ejected crust. 

We find that although a torus and shell may seem quite sim- 
ilar, the effects of geometry on bursting and period glitches is 
significant. As such, in this paper, we attempt to describe the 
fracturing of the torus due to shear forces caused by differential 
rotation within the torus. Moreover, we show how this leads to 
the accretion of the inner edges of the torus as it is slowly per- 
meated by the star's magnetic field, leading to X-ray bursting. 
The nature of this bursting, as seen in the AXP 1E2259+586, 
can be explained within the framework of our model by appeal- 
ing to a degenerate torus with a high metallicity, which exhibits 
changes in it's ionization fraction as it is bombarded with X- 
rays during bursts. It may be the case this emission, due to the 
presence of a torus, is insignificant compared to emission from 
vortex expulsion. However, we feel that in at least the cases of 
the two oldest AXPs (1E2259-H586 and4U0142-H615) the torus 
emission can exceed emission from vortex expulsion. 

This paper is presented as follows: In §|2]we review the 
transport of angular momentum leading to the torus forma- 
tion and it's evolution. In § |3] we discuss the quiescent phase 
due to accretion from the non-degenerate atmosphere of the 
torus, followed by the bursting phase in § |4]due to the star's 
magnetic field penetrating the torus's inner edge. Following 
that, we apply our model in §|6]to the AXPs 1E2259-I-586 and 
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4U0142+615, and in §|2]we provide a brief discussion of the 
implications and predictions of our model. 



2. Propeller regime and torus formation 

As described in Keranen et al. 2005, during the quark-nova the 
core of the parent object converts to (u,d,s) quark matter and 
becomes suddenly compact. As the core contracts it becomes 
detached from the outer crust. Consequently the outer crust is 
left at a radius larger than the surface of the newly born quark 
star The fate of this crust depends on the amount of energy 
released during the quark-nova, and its coupling to the crust. 

The metal-rich ejecta from the Quark-Nova (QN) is de- 
generate (described in OLNI), and remains so as it expands 
out to the magnetic equilibrium radius, R^, which is the ra- 
dius where the ejecta's gravitational pressure balances the star's 
magnetic pressure. Once there it will have the form of a shell, 
and will have expanded only to a thickness of, A/?ni/^m = 
1.2 X I0^^m[^*, where the shell's mass in terms of IO^^Mq is 
m_7. Provided that the quark star's (QS) magnetic field is strong 
enough to support this shell, then it will cool rapidly remaining 
degenerate, as shown in OLNI. 

If this magnetic equilibrium radius is larger than the coro- 
tation radius (which should be possible in at least some in- 
stances given the uncertainty in initial conditions of the quark- 
nova compact remnant) the propeller mechanism will take ef- 
fect (Schwartzman 1970; Illarionov & Sunyaev 1975), deflect- 
ing the degenerate shell into a torus on the equatorial plane, 
which is the situation considered in this paper 

Using an angular momentum conservation argument, we 
can estimate the location of such a torus by writing 



(1) 



where QQs,i is the quark star's initial period, Ri is the equa- 
torial location of the torus in a Keplerian orbit, and Qk = 



4 



GMqs/R^ is the Kepler rotation frequency. Since all of the 
shell material is propelled away at the equator with constant 
specific angular momentum leading to a thin torus. Applying 
the equation for the magnetic radius (Eq. 1 in OLNI), implies a 
torus radius of. 



TJt ^ 15 km 



d4 nil 
°0,15"qS,10 

" -V'^QS, 1.4^ i,mi 



(2) 



where the QS's birth period, Pi,ms, is given in units of millisec- 
onds, while the QS's surface magnetic field strength at birth, 
So,i5, is in units of 10'^ G. The radius and mass are ^qs,io 
in units of 10 km, and Mqs,i.4 in units of IAMq respectively. 
However, in order to form a torus we require enough angular 
momentum transfer to guarantee Rt > 7?nj. This translates into 
an upper limit on the initial period of. 



Pims < Plim ~ 2.5 ms 



r3/2 „9/2 



QS.IO 



■'"qS,!.. 



(3) 




0.001 



Fig. 1. X-ray luminosity is plotted against period derivative ac- 
cording to Eq. (30) in OLNI. The upper and lower solid lines 
represent a magnetic to X-ray conversion efficiency, 77X, of 1 
and 0. 1 respectively, and the dashed line is the luminosity av- 
eraged over all viewing angles. AXP 1E2259+586 and AXP 
4U0142-I-615 have luminosities higher than what is predicted 
in OLNI from magnetic field decay alone. The higher luminos- 
ity we attribute to steady accretion from a torus. 



Moreover, the duration of the propeller phase can be ap- 
proximated to be fpi-op ^ Rt/vprop. With Vpi-op ^ Oqs x R^ this 
gives 



1 



1 



tprop. 



0.15 ms 



i,ms t.l5 



'QS,i 



M, 



1/4 



(4) 



QS,1.4 



2.1. Torus fracture from tidal shear 

After vertical expansion the torus radial extent (eq. lA.6l ) is still 
small. However, the Keplerian velocity within the torus varies 
from vk,i - -\JGMqs /Rt at the equator to vk,2 - -\/GMQs/1.7Rt 
at z = ±Rt (see eq. IA.3b . resulting in a large shear. This shear 
force acts only in the azimuthal direction, or in other words the 
shear operates only between surfaces of differing Vcirc. Thus, 
the main part of the torus where density is highest will fracture 
into subsequently large cylinders, or "walls" of thickness Ar„ 
and mass m„. Hydrostatic equilibrium along the surfaces of the 
walls guarantees that each wall will remain in high pressure 
contact with the next. In other words we do not expect a dis- 
continuous set of walls but a rather walls separated by melted 
(non-crystalline) material, or basically a degenerate fluid. 

To determine the thickness of these walls, one appeals to 
the ultimate strength of iron, which is f/ = i/'£'b, where Eb ~ I 
MeV is the binding energy per nucleon of the torus material, 
and t// is the strain factor t// ~ 10"^ for Iron (e.g. Halliday 



Ouyed et al.: Tori in Quark-Nova remnants: application to AXPs 



3 



& Resnick §13.6). By equating the work performed by the 
shear forces W = F x Ar^, (where the tidal force per unit 
area is F = -IGMQ^Ar^/ r^,) with the torus binding energy, 
£B,Mevwiw/(56mH), One finds a typical wall thickness to be 

Ar„ ~ 400 cm fperlj^- Here fy, = '/'is^^'aMev/^Q's.i^' 
rw,i5 is the distance of the wall from the center of the quark star 
in units of 15 km. Thus, as a fiducial example, the innermost 
wall (r„ = Ri) has a thickness of. 



Ar„ ~400cm/Fe/?, 



3/2 



(5) 



The number of walls in the entire torus is then dependent 
on the thickness of the torus, and can be found to be. 



A^w ~ 2000- 



p5/4 
- keV 



3/2 
t,15 



(6) 



Also, the mass of the innermost wall in terms of the torus equi- 
librium temperature (equation lB.13b is. 



15/8/. nl/'t 

'0.1 yrs 



1/4 



mw ^ 6.5 X 10""M, 



(7) 



where /^atm is the mean molecular weight of the torus atmo- 
sphere discussed in Appendix A and in the following section. 
The time dependency in the equation above is due to the aver- 
age density of the torus decreasing in time as the torus spreads 
radially in a manner given by Eq. ( |A.6t . There is also a small 
extra dependency in time due to the torus mass gradually de- 
creasing from torus 's atmosphere being accreted during the 
quiescent phase, thus reducing the total torus mass as given 
by equation ( fTOb . 

3. The quiescent phase 

As explained in Niebergal et al. (2005), the X-ray luminosity 
during the quiescent phase in our model is due to vortex expul- 
sion. The magnetic field contained within the vortices is also 
expelled, and the subsequent magnetic reconnection leads to 
the production of X-rays with a luminosity as given in OLNI, 

Lx ^ 2.01 X 10^^ erg &-^t]^P\^ , 



(8) 



where rjx is the efficiency parameter inherent in the coversion 
from magnetic energy to radiation. Figure 1 (reproduced from 
OLNI) shows our model of luminosty evolution as compared to 
the quiescent X-ray luminosities of AXPs, SGRs and XDINSs. 
It is clear that AXP 1E2259 and AXP 4U0142 show an excess 
luminosity as compared to SGRs/AXPs born with a shell. This 
we attribute to continuous accretion from the torus as explained 
below. 

3.1. Steady accretion from the torus 

At the edges of the torus, using equation ( IA.4l i. the den- 
sity is low enough such that the torus material becomes non- 
degenerate, which creates a thin metal-rich atmosphere. A frac- 
tion of this atmosphere is constantly accreted onto the QS, 
where it is converted into CFL quark matter, releasing the ex- 
cess energy as radiation. This radiation in turn heats up the 
torus, thus altering the size of the non-degenerate atmosphere. 



For a given temperature, material from the non-degenerate 
torus atmosphere slowly leaks out at a rate, ih. For any non- 
zero ih, the accretion luminosity heats up the torus, which also 
cools rapidly as a blackbody. Therefore, one can determine the 
equilibrium temperature given by Eq. ( IB.13I ), 



Teq ~ 0.52 keV 



^atm^QS,1.4 



(9) 



This also implies a continuous or equilibrium accretion rate 
from the torus atmosphere of. 



2.3 X 10^'' gms 



A*atm-'"QS,I.4 

Thus, the corresponding equilibrium accretion luminosity is 

_1 %.l^t,15 



(10) 



Lace - J]meqC^ ~ 2.2 X lO"'^ erg s 



(11) 



QS,1.4 



This defines the quiescent phase luminosity in our model 
for a QS with a torus. From this, a blackbody temperature from 
the accretion hot spot on the star's surface can be found to be. 



rBB,eq~0.92keV 



d3/2 



/^atm^BB,5^QS,1.4 



(12) 



where Rbb.s is the accretion region (i.e. the polar cap since most 
of the accreted material is channeled towards the poles) in units 
of 5 km. 

3.2. The effect of mean molecular weight, ju 

Ionization of the torus atmosphere is due to X-ray photons 
from accretion onto the poles as mentioned above. Iron ion- 
ization equilibrium has been discussed at length in Arnaud 
& Raymond (1992), who show (see their figure 10) the most 
abundant charge states versus temperature of iron. For the equi- 
librium temperatures in the torus (Eq. |9]l, we estimate that the 
charge state is situated between the L and M shells, which cor- 
responds to yUpe values of ~ 3.3. For comparison fully ionized 
iron yields //pe ~ 2.07. 

Therefore during the quiescent phase (p. ~ 3.3), our model 
predicts the following parameters and observables. 



1.9 X 10" gms 



3 /j6 

13 -I -'O.r^t.is 



■"^05,1.4 



(13) 



~ 0.085 keV ' ' 

Mqs,1.4 

d3/2 

rBB,eq ~ 0.15 keV 

T 1 T irv33 -1 'o.rS.is 

Lace ~ 1.7 X 10 " erg s — . 

■^^08,1.4 

Clearly, in the case where a torus is formed around a QS, 
as opposed to the shell case presented in OLNI, the quies- 
cent phase is dominated by continuous accretion luminosity, 
not by vortex expulsion. The observed quiescent luminosities 
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of AXP 4U0142 3.3 x lO^'* ergs-') and AXP 1E2259 
3.8 X 10^"* erg s"' ; see Table I in OLNI) implies the location 

of these tori to be at roughly Rt ~ 25 km. These two candidates 

are discussed in more details in § |6] 

In principle, for a given source (i.e. fixed Rt) a change by 

a factor of a few in the mean molecular weight of the torus 

atmosphere implies large changes in the accretion luminosity 

from the torus since 



t,eq 



' BB.eq 



Lis 



-3/2 
^atm. 

-3/2 
/^atm. 

/^atm. ■ 



(14) 



For example if the degenerate atmosphere is suddenly heated as 
to become fully ionized during a bursting episode (discussed in 
the next section), then yUatm would drop from 3.3 to 2. 1 yielding. 



eq 

T^BB.b 



2.0 



2.0 



(15) 



' BB.eq 



^b,a( 



16.4. 



where the subscript "b" stands for bursting. During the bursting 
phase then, the torus equilibrium temperature and the black- 
body temperature of the hot spot on the star would have in- 
creased by a factor ~ 2 while the torus accretion luminosity 
would increase by over an order of magnitude from its quies- 
cent value. The mean molecular weight can also change due to 
a change in composition leading to enhanced temperature and 
luminosity (see ^4.41 1. 



where, Ar^, is the wall's thickness (see Eq.|5]l, and equation [A.8l 
was used. As the torus spreads in time its density decreases as 
f-i/2 (fi-om Eq. lA.Sb . thus weakly decreasing the penetration 
timescale as the source ages by, tb oc r'^'^. 

4.2. Accretion of the wall 

Once the magnetic field penetrates the innermost wall, the 
poloidal component of the dipole field is wound up introduc- 
ing an azimuthal component B^. The amplification of the mag- 
netic field can be found from 5B/5f = V x (v x B) which yields 
dB^/dt ~ QkBp, where Qk = y/(GMQs/Ri) is the Kepler fre- 
quency and Bp is the poloidal component of the QS's dipole 
field. That is, the magnetic field within the wall will be am- 
plified following the relation B^ oc Q^Bpf, until the magnetic 
energy is comparable to the kinetic energy of the wall, at which 
point it dominates the dynamics. This rapid build up of B^ oc- 
curs roughly on the timescale of milliseconds. 

The resulting magnetic torque (oc R^B^Bp) spins-down, the 
wall and spins-up the star. Or in other words, the magnetic 
torque is essentially separating the wall from the main body of 
the torus and transferring it's angular momentum to the star, in- 
ducing a rotation period anomaly (called a "glitch"; discussed 
further in §|5]i. This causes the wall to eventually co-rotate with 
the QS. 

The wall's initial mass at the time of detachment from the 
torus, m^ i, is given by equation (|7]i (using /iatm = 3.3), 



6.1 X 10" "'Mp 



ff^m^iR^ 

5/4J/2 



1/4 



(18) 



yrs 



4. The bursting phase 

In our model, the bursting phase is initiated when the innermost 
wall of the torus is magnetically permeated such that it detaches 
from the torus itself, where it is then accreted rapidly onto the 
QS. 



4.1. Magnetic Held penetration 



Recall from § 12.11 that the torus is composed of subsequently 
larger walls, made of a solid degenerate matter, and separated 
by a degenerate fluid. As such, the torus's inner wall will be 
slowly penetrated by the QS's magnetic field on timescales de- 
termined by the induction equation. 



dB 

dt 4-ncr 



(16) 



Here, cr = ne,the^/le/('«eVims) and Ag = 1 /(ne,thO"T), where ng.th 
is the number density of thermal electrons in the degenerate 
matter, crj is the Thompson scattering cross-section, and the 
root-mean-square electron velocity is vvms = Cj. Therefore, the 
time needed for the magnetic field to penetrate to a depth of 
Arw into the torus is. 



Tb ~ 82 yrs 



Fe"t,I5 



(17) 



Thus, as the the gravitational pressure of the wall is much 
weaker than the magnetic pressure from the QS, the co-rotating 
wall is trapped between the torus and the pressure of the un- 
derlying magnetic field. Furthermore, the effective gravity act- 
ing on the wall's edges also guarantees that the wall is pres- 
sure confined vertically. Moreover, by equating the magnetic 
pressure to the wall's pressure, B^/(87r) = KPw'^^^ where k = 
1.24 X lO' Ve^^^' we find the density of the wall after detach- 
ment to be. 



1.38 X 10** gm cm" 



^3/2^9/2 



QS.IO 



Ri: 



,9/2 
t,15 



(19) 



where the subscript "w" implies a co-rotating wall and depicts 
the fact that the detached wall is spun-down to co-rotate with 
the star. The value of the magnetic field in equation above is 
evaluated at detachment and is given in our model as B^ - 
y/3KBPP where kb = 8.8 x 10^** G- s"' (see Niebergal et al. 
2006; and also OLNI). 

Given this situation, the wall will lose mass by thermal 
evaporation of it's non-degenerate atmosphere, thus getting 
thinner in time while keeping the same density and pressure 
(i.e. Ar„ - 2KRim^(t)l p^). By appealing to equations dC.ll l 
& (|C.2t in the appendix we arrive at the following equation 
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for the decrease in mass of the wall (expressed here in units of 
IO^'^Mq) over time, 



dm 



w,-10 



-7.7 X 10" 



-dt. 



Thus, the wall is accreted in time at a rate of. 



WIw,i , f 

= 1 + — , 

Tw 



(20) 



(21) 



where a characteristic time for wall consumption is defined to 
be, 



Tw =0.4 hours 



/^w,atmfis_i4^QS,io / 10"'°M( 



'70.1^; 



.5 

t,15 



(22) 



The expression for t„ implies that the more massive the wall 
at detachment the faster it gets consumed. Equation ( |2TI ) shows 
that 99% of the wall is consumed in 100t„. The accretion of the 
wall gives a much higher accretion rate (see eq.(B 1)) compared 
to the quiescent phase accretion from the torus atmosphere. The 
resulting equilibrium temperature (eq.(B2)) is high enough to 
fully ionize the wall atmosphere, resulting in a yUw.atm ~ 2.1 
while the wall is being consumed. 

4.3. The first burst component: the steep decay 

The first component in our model is defined by consumption of 
the detached wall. Using equation (ISTT l. the corresponding lu- 
minosity due to accretion of the wall from its atmosphere varies 
in time as. 



Lw (t) - Tjm^c 



(23) 



(1 + i)' ' 

where the initial luminosity due to the accretion of the wall is 
r]o.\ ( '«w,i 



^^,.5X10 ergs — (^^) , 



(24) 



and Tw,hrs, is expressed in hours. The wall temperature as given 
by equation (B2) starts at a few keV and drops as mw(f)'^^ until 
it becomes negligible compared to the steady emission due to 
the torus accretion. See Figure Q for a comparison with ob- 
servations. 

During the entire accretion process the amount of energy 
released will be. 



-f 

Jo 



Lw (0 dt — tLw.o 



(25) 



5x lO^"* erg 770.1 



lO-i^Mf 



4.4. The second burst component: the slow decay 

During the burst, the torus is irradiated by MeV neutrinos and 
photons causing the dissociation of the iron in the torus. For a 
typical total energy of about 10'''' erg (~ 10^^ MeV) released 
during the consumption of the detached wall (see Eq. |25j, we 
estimate about 1 iron nucleus dissociated per MeV, or a total 



of ~ 10"*^ dissociations. Hydrostatic buoyancy then causes the 
nuclei with Z~ 26/2 =13 to float to the top of the torus onto 
the non-degenerate atmosphere. Following the wall evapora- 
tion and torus irradiation, the torus atmosphere is now mostly 
composed of nuclei with Z ~ 13 instead of Z = 26. This result 
in decreased //atm since the lighter nuclei are fully ionized (they 
have lower atomic energy levels by a factor of Z^ ~ 4) causing 
yUatm to drop as low as fintm.b ~ 2.1 (in the extreme case of a 
pure light nuclei atmosphere) from //atm.q ~ 3.3. 

Over time, /iatm increases again to 3.3 as the 10"^^ light nu- 
clei in the atmosphere are slowly depleted by accretion while 
the atmosphere gradually becomes enriched to iron again. 
A simple model for light nuclei depletion yields dXn/dt - 
-Xo/to where X13 is the fractional abundance by mass of 
the light nuclei in the atmosphere and T13 is the exponential 
decay timescale. This timescale is the mass in light nuclei, 
OT13 ~ 3 X IO'^^oth, in the atmosphere divided by the depletion 
rate given by tht, or. 



'«13 , , 
= — ~ 13 year- 

'?o.i^t. 



^QS,1.4 



(26) 



found using an average of 2.5 for the mean molecular weight 
per electron. 

We model the mean molecular weight using two species: 
partially ionized iron (A - 56, Z = 26, Zg - 16), and fully ion- 
ized light nuclei (A = 28, Z = Zg = 13) where Zg is the number 
of free electrons per nucleus of charge Z. The corresponding 
molecular weight of the mixture is /i"' = + Zej)/Ai 

which after substitution for Xuit) and X26 = 1 - ^13 yields 

1/1 1 . ^27) 



A'atm(0 /^atm,q \/^atm,b /^atm.q ^ 

With the use of equation (fTTT i. the luminosity due to accretion 
during a burst is then. 



ib,a 



-^q.acc 



l+( 



^atn 



A'atm.b 



l)e 



-tiTu 



(28) 



where Lq ^cc is the accretion luminosity during the quiescent 
phase when the torus mean molecular weight is yUatm.q ~ 3.3 
(see Eq. [TTI ). A comparison with observations using this is 
shown in Figure (|2]l. 

5. Changes in Rotation Period 

5.1. Instantaneous spin-up during a burst 

After the star's magnetic field has permeated the innermost wall 
of the torus, the wall will detach and begin to co-rotate with the 
star During this transition, by conservation of angular momen- 
tum, the star's period will decrease by an amount given by the 
expression, 

AP MM 5 m„ Q., 

which can be instead expressed in terms of the star's rotation 
period (in units of 10 seconds) by, 

— -3.8x 10-^^^^^ I TTTTirrri- 1 ■ (30) 



(29) 



Mqs,i.4«2^_jo\10-i"M, 
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It is worth pointing out that the period glitch is not neces- 
sarily simultaneous with the burst, since there may be a period 
of pressure adjustment before the wall begins being consumed. 
Thus, a clear prediction is that QSs born with a torus instead of 
a shell will show glitch activity just prior to their bursts. 

5.2. Persistent spin-down following bursting 

A portion of the non-degenerate wall atmosphere matter is kept 
in co-rotation by the magnetic field out to the light cylinder, 
which provides an efficient mechanism for removing angular 
momentum from the system, via the propeller The angular mo- 
mentum per unit mass lost at the light cylinder is c^/Qqs, which 
enhances the spin-down rate of the quark star during bursts to, 

P_i3,b - 2m„,,io^'io ■ (31) 

Here, the wall's evaporation or mass-loss rate, niH.jQ, is given in 
units of 10'" g s and the enhanced spin-down rate is in units 

of 10-1^ s s-i. 



Using the observed quiescent phase luminosity of AXP 
4U0142+61 (Lace = 3.3 x lO^'* erg s"'), we find ^ 25 km. 
From this, our model yields the equilibrium temperature of the 
torus, blackbody temperature of the emitted X-ray luminosity, 
time in years between bursts, and the age of the system: 

Rt - 25 km 

Teq ^ 0.14 keV 

TBB.eq - 0.42 keV 

tb - 150 years (32) 

fage ^ 4.7 X 10^ years (33) 

A 5 km blackbody emitting radius was assumed to get the 
blackbody temperature. When a burst occurs in the future, 
equation i25[ gives us the wall mass then equation (l30l l yields 
the expected glitch, AP/P, in our model. Using R^ ^ 25 
km we find AP/P ~ -9 x lO"^(mw,i/lO"'°M0). The torus 
mass can then be derived from equation (fTsT l to be nit - 
2 X lO-''Mo(OTw,i/lO"'"M0) which is of the order of IQ-^Mq 
for our fiducial wall mass. 



6. SGRs and AXPs in our model: Case study 

Here we specifically focus on two AXPs, namely 1E2259-I-586 
and 4U0142 -1-615, (see Table 1 in OLNI for more on their ob- 
served features). These two candidates, as can be seen in Figure 
1, show luminosity during their quiescent phases above the lu- 
minosity predicted by vortex expulsion. In our model this ex- 
cess luminosity is due to the presence of a torus, which is con- 
stantly being accreted. Evidence for the presence of a torus 
has indeed already been detected around AXP 4U0142+615 
(HuUeman et al. 2000; Wang et al. 2006). 

6.1. AXP4U0142+615 

Observations of pulsed optical emission from AXP 
4U0I42H-61 has led to the application of the irradiated 
disk model, wherein optical and infrared luminosities 
are reprocessed emission from an X-ray irradiated disk 
(i.e. Hulleman et al. 2000a). Most irradiated disk models 
predict an optical emission larger than what is observed in 
4U0142H-61 (Hulleman et al. 2004), so to compensate one 
either assumes a larger inner disk radius or a smaller outer 
radius. The former implies the absence of a hot inner region 
which is inconsistent with observations, leaving the option of 
an outer disk radius very close to the inner radius. Thus, one 
has the possibility of a very thin passive disk, which we argue 
fits all of the features of the degenerate torus as described in 
this paper. Normally a very thin disk is considered unrealistic 
due to heating by viscous dissipation (Hulleman et al. 2000b), 
however, in our model viscous dissipation is negligible due to 
degeneracy in the disk. 

Moreover, if the blackbody component is to be believed, 
then the emission region is confined to a radius of 12 km (White 
et al. 1996; assuming a distance of 5 kpc), ruling out some mag- 
netospheric emission models. In our model, an emission region 
of this size is predicted, and confines parameter values that fit 
well with other observations (see Eq.fT2ll. 



6.2. AXP1E2259+586 

During the June 2002 outburst (Woods et al. 2004), AXP 
1E2259H-586 displayed over 80 X-ray bursts for approximately 
4 hours with bursts ranging in duration from 2 ms to 3 s. We re- 
fer to this as the first component which decayed approximately 
as a power law in time, oc r"^-^. Enhanced flux was observed 
over the following year, referred to as the second component, 
which decayed according to a power law in time oc f""-^^. The 
X-ray properties of the bursts are very similar to those seen in 
Soft Gamma Repeaters (SGRs). This is natural in our model 
since SGRs and AXPs involve the same underlying engine, 
namely a quark star accreting from a degenerate shell/torus. 
We thus propose that the only difference between SGRs and 
the two AXPs mentioned in this paper, is that the SGRs were 
born as QSs with a slower rotation period, therefore possessing 
a co-rotating shell (as discussed in OLNI), rather than a thin 
torus. 

Observations give the quiescent phase blackbody tempera- 
ture of about 0.42 keV. Given the observed quiescent phase lu- 
minosity our model implies Ri ^ 25 km. From this, our model 
yields the equilibrium temperature of the torus, blackbody tem- 
perature of the emitted X-ray luminosity, time in years between 
bursts, and the age of the system: 

Rt ^ 25 km 
Teq ^ 0.14keV 
T^BB.eq - 0.42 keV 

Tb ^ 100 years (34) 
fage - 1.5 X 10^ years (35) 

During the bursting phase and at the onset of the outburst 
the measured temperature increased to 1 .7 keV, before decay- 
ing to 0.5 keV within the first few days. This is associated 
with the wall accretion episode in our model. Gavriil et al. 
(2002) reported that they simultaneously observed increases 
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Fig. 2. Shown is the two components of the burst history for 
AXP 1E2259, where the soHd line is the steep decay (Eq.l23Tl. 
the dash-dotted line is the slow decay (Eq. |28] ). the dotted line 
is the combination of both the steep and slow decay, and the 
dashed line is the quiescent phase luminosity (Eq.fTTTl. The data 
points are taken from XTE and XMM observations (i.e. Woods 
et al. 2004). 

of the pulsed and persistent X-ray emission by over an order 
of magnitude relative to quiescent levels. This increase is also 
consistent with our model because, as discussed in § 13.21 re- 
duction in the mean molecular weight of the torus atmosphere 
during the burst causes a significant increase in the temperature 
(Tbb ^T^^^ ) and luminosity (i.e. L cc yu"^). 

Further observations of AXP 1E2259+586 show that it un- 
derwent a sudden spin-up (AP/P = -4 x 10 ''), followed by a 
factor of 2 increase in spin-down rate, which persisted for more 
than 18 days. The spin-up is consistent with the torque due to 
wall detachment from the torus (see ^S.ll l. Using equation dSOl l 
and the measured period glitch, one gets for the mass of the de- 
tached wall, mw - 0.5 x 10 '"M©, which is in good agreement 
with the expected wall mass based on shear forces on the torus 
(Eq. [TST l. Moreover, using this wall mass, our model predicts 
(using Eq.l25Tl the energy in the burst to be 2 x 10^^ erg, which 
is an order of magnitude greater than the observed burst, im- 
plying that observations likely missed the early phases of the 
burst. In addition the age and the wall mass allow us to calcu- 
late the torus mass from equation (fTST l to be of the ordei^ of 
nit ~ 2.2 X IO^^Mq. 



^ Interestingly, Wang et al. (2006) estimate the mass of the disk sur- 
rounding 4U0142-I-615 to be within the same range as the torus mass 
in our model. We are tempted to speculate that the discovered debris 
disk is in fact the degenerate torus described in our model. Planet for- 
mation around such a metal-rich degenerate torus has been suggested 
and discussed in Keranen&Ouyed (2003). 



The characteristic timescale for the two components in our 
model in the case of AXP 2259-1-586 are given by equations 
(Hall and (ESJ, 

Tw ^ 0.11 hours 

Ti3 ^ 141 days. (36) 

The time needed to consume 99% of the wall is ~ IOOtw ~ 
10 hrs, consistent with the observed June 2002 outburst. 
We note that the observed enhanced spin-down of AXP 
1E2259+586 lasted for ~ 18 + 6 days. Naively, we expect the 
enhanced spin-down episode to coincide with the wall accre- 
tion event. However, the propeller torque could in principle act 
for longer period if there is more matter supplied beyond the 
wall accretion episode. Shown in Figure (|2]l are the two pre- 
dicted decay components as compared to the observed outburst 
of AXP 2259-1-586. The two components in our model are given 
by equations ( |23] l and ( |28] l. 

The second component in our model is consistent with the 
data, if we adopt yUatm.b ~ 2.8, as shown in Figure |2] Adopting 
yUatm.b ~ 2.8 suggests that the composition of the torus atmo- 
sphere following irradiation is slightly different from a pure 
light nuclei atmosphere. The predicted decay of the first com- 
ponent is slower than the observed one. We argue that a faster 
decay rate could result if one were to include "squeezing" pres- 
sure exerted on the wall by the dipole magnetic field. Such an 
elaborate description of the wall consumption could also lead 
to sporadic episodes of high accretion which could account for 
the ~ 80 X-ray bursts observed in this source. A minimum 
timescale for such accretion episodes, corresponding to wall 
material free-falling onto the star along the magnetic field lines, 
is 

face = — ~ 1 ms — ^ 2.15 ms , (37) 

™QS,I.4 

which matches the lower limit in the 2 ms to 3 s range in bursts 
duration observed; Vif - y/lGM/Rt is the free-fall velocity. 

7. Discussion and Conclusion 

With this paper and OLNI, we have completed our investiga- 
tion of the fate of the ejected crust during a quark-nova event. 
We conclude that, depending on the unknown initial conditions 
of the quark-nova compact remnant (specifically the initial spin 
period, magnetic field strength, and ejecta mass), the ejected 
crust (from the parent neutron star) will form either a shell or 
torus. 

Given that the initial spin period that separates torus for- 
mation from shell formation is close to values inferred the- 
oretically (Staff et al. 2006), one should expect that some 
SGRs/AXPs possess a shell and some a torus. In OLNI we 
concluded that most of the SGRs/AXPs are best modeled by 
appealing to the shell picture. In this paper, we have shown ev- 
idence that two AXPs, 4U0142+61 and 1E2259+586, are bet- 
ter understood in the torus picture. These AXPs show X-ray 
luminosity in excess compared to what is expected from vor- 
tex expulsion, given their rotation period and period derivative 
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(see Eq.[8]). We argue that these sources were merely born with 
initial periods below the limiting period given in equation (O. 

Most of the dynamics and all of the emission mecha- 
nisms applicable to a torus should also apply to the shell case. 
However, in most SGRs/AXPs, all emission is dominated by 
flux expulsion, so these effects are unimportant except in AXPs 
1E2259-I-586 and 4U0142-t-61. Thus, it may be the case that 
other SGRs/AXPs also possess a torus. 

The recent work by Durant & van Kerkwijk (2006) in- 
dicates that all AXP luminosities are roughly equal to 1.3 x 
lO^^ergs"'. Although these improved luminosity estimates 
may indicate a constant luminosity for all AXPs (and pos- 
sibly the same emission mechanism), these improved lumi- 
nosities actually better agree with our theoretical fit (c.f. dot- 
ted line in Fig. 1). Regardless, one still needs to consider the 
reason for the lower period derivatives of AXPs 1E2259-I-586 
and 4U0142+61. Equation ( fTTT ) predicts that the quiescent ac- 
cretion luminosity in our model is independent of the period 
derivative, so when this constant accretion luminosity becomes 
roughly the same as the spin-down luminosity (eq.[8]l, we can- 
not distinguish whether the object has a shell or torus. 

Although accretion and propeller models have been con- 
sidered in the past (eg. Mereghetti & Stella 1995, Mosquera 
Cuesta et al. 1998, Rothschild et al. 2002), as well as models 
involving exotic compact stars (eg. Usov 2001), our model is 
the first to use both in the context of a CFL quark star to explain 
both bursting and evolution. Also, there are clear advantages to 
our model; i) whereas normal accretion disks are generally not 
able to survive the intense gamma bursts, the dense ejecta in 
our model can not only survive but also reacts in a manner that 
provides a good fit to burst profiles. Also, ii) the presence of a 
CFL quark star has the advantage of being able explain both pe- 
riod clustering and long-term spin-down by appealing to vortex 
expulsion, as well as possessing the environment necessary for 
the type of disk we describe, which accounts for bursting be- 
havior; and iii) the CFL quark star is extremely volatile, which 
makes achieving the large burst energies does not require ex- 
cessively large magnetic field strengths. 

We realize that although it may appear elaborate, the wall- 
Uke structure is a result of the first-order treatment of the inter- 
play between gravity and differential rotation. In reality we ex- 
pect a more complex formation and evolution of the torus, but 
we feel the simplified geometry (i.e. "walls") helps understand 
the essentials. Furthermore, we have not studied the stability of 
this equilibrium structure. We feel that stability could be an is- 
sue, however, the model's ability to account for observations of 
AXPs indicates the structure may indeed be stable and similar 
to what we describe. 

A further impUcation of our model is that the initial spin 
period of the quark star is determined by observational param- 
eters. Equations (2) and (3) in Niebergal et al. (2006) imply a 
conserved quantity in our model given by 



PB^ = PoBf. 



(38) 



Equation above is valid if accretion and propeller torques are 
small compared to magnetic spin-down torques averaged over 



long timescales. Combining this with our equation (|2| yields 
the quark star birth period 



5 ms 



1/2 
-13 



dI/4 1/2 
^t,15'«-7 



(39) 



For example in the case of AXP 2259-H586 for which we have 
derived the torus mass and radius we find a corresponding birth 
period of ~ 5 ms for a torus mass of m_7 ~ 2. This is below 
the limiting period, given by equation (O, implying torus for- 
mation, Piim ~ 8 ms, providing a self-consistency check of our 
model. Additionally, from the observational point of view, the 
lack of energetic SNR surrounding some AXPs and SGRs such 
as AXP 2259+586 (Vink&Kuiper 2006) impUes a birth period 
greater than about 5 ms. 

It is not clear if the passive disk discovered around AXP 
4U0142 is the torus we describe here or just a fall-back disk, 
as in the case of rapidly rotating progenitors, which could have 
formed directly from the collapsing SN envelope (before the 
QN explosion). However, the low optical emission of the dis- 
covered disk, is highly suggestive of a degenerate torus as we 
described in our model. Further observations of SGRs/AXPs 
should be able to distinguish between the disk and the degen- 
erate torus, and would provide an excellent test for our model. 

Acknowledgements. This research is supported by grants from the 
Natural Science and Engineering Research Council of Canada 
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Appendix A: Torus geometry and evolution 

The formation of the torus involves a vertical expansion at the speed 
of sound until hydrostatic equilibrium is reached. For a completely 
degenerate and relativistic gas the sound speed is (see OLNI), 



= 1.06 X lO'VIo^ cm^/s^ 



(A.1) 



where pio is the torus density in units of 10'" g/cc. 

While the material expands vertically, momentum conservation 
(i.e. Vii(Rt)Rt = Vci,c.(z, r)r, where vk is the Keplerian velocity) causes 
the material to move radially at the same time. Combined with the 
circular orbital motion at height z. 



GM, 



QS 



cos(ff) , 



(A.2) 



where cos(£') = r/(r^ + z^)''^ and (r, z) are cylindrical coordinates, the 
resulting shape of the expanding torus is given by, 



(A.3) 



From equation ( |A.3! ) above one can see that the torus flares into a con- 
ical shape. Integrating the equation of hydrostatic equilibrium along 
the flared surface of the torus from {z = Q,r - Ri) to rsph. = Vr^ -i- z^ 
gives the density profile 



pinzY'^ =p(Rt,z = 0) 



1/3 , 



GM, 



QS 



4k 



Jo ~- 



+ z' 



ds , 



(A.4) 



where (f> is the angle between the tangent to the surface (given by 
Eq. IA.3t and the radial direction {ds is the path length element along 
the surface). 
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Using equation JA.3t , the resulting structure extends from (r, z) ~ 
(lARi,-Rt) to (r, z) ~ {l.4Rt,+Rt) and is independent of the torus 
mass. In simple terms, the effective gravity (g^ff = -GMgs/Cr^ + z^)) 
is small from z = to z ~ Rt while above where the surface 
is oriented nearly radially, the effective gravity is strong resulting in 
a rapid decline in density. The timescale for vertical expansion is 
Rt/Cfi ~ (1-5 ms) X i?t,i5/p,'io where p, lo is the torus density in units of 
10'° g cm-\. 

The radial thickness of the torus, Ar,, is governed by the angu- 
lar momentum transfer due to viscosity. The viscosity due to particle 
collisions in a degenerate gas is estimated using that for an ideal gas 
(e.g. Lang pg. 266fl 

v = 6x 10"r^'^'vcm's-' . (A.5) 

Equation l lA.St shows the strong dependence of the viscosity on the 
temperature, thus on cooling. The radial expansion of the torus can 
be derived using d ((Ar)-^ /dt = v. Due to the short timescale associ- 
ated with the propeller phase, (eq.|4j, the radial width prior to vertical 
expansion is 



(^'')prop (^^prop) 



7^5/4 „3/8 
1-7 1 /-v4 McV"l,15 

1.7 X 10^ cm- 



^"qS,1.4 



(A.6) 



As shown in OLNI, the shell temperature at formation is of the order 
of 1 MeV. Blackbody cooling of the shell during the propeller phase 
yields a temperature of the propelled matter decreasing with time from 
an initial ~ 1 MeV to a final ~ 60 keV. Including cooling during 
the propeller phase slightly affects (Ar)prop ; the constant in front of 
equation l lA.6t reduced to 1.2 x 10"* cm. This is understandable since 
the torus spreading is dominated by the early hot part of the propeller 
phase. During the vertical expansion, with an initial temperature of 
~ 60 keV, viscous spreading results in an additional (Ar)exp ~ 10^^ cm. 

Following the vertical expansion, accretion from the edges of the 
torus begins. During this time the torus is also cooling as a blackbody, 
and quickly reaches an equilibrium temperature, T^q, in the keV range 
as given in Appendix A (ea.l lB.13} ). Subsequent viscous spreading 
is then determined by T^q. Thus the spreading for some equilibrium 
temperature is. 



r5/4 1/2 
' keV'yrs ' 



(A.7) 

where fyrs is the age of the system in yearfl and Ti^ev is the equilibrium 
temperature of the torus in keV. The torus average density decreases 
in time as. 



47rR^Art 



: 9.1 X 10^ gcm"^' — 



m_7 



pi 7^5/4,1/2 



(A.8) 



In general at a given temperature the maximum density, Pnd.q^ of 
the torus below which the matter is non-degenerate is found by setting 
Tt = rpcrmi - 123.6 MeV/j(r"''pjg' (see appendix A), or, 



p„d - 230ytie gem ^tI^^ . 



(A.9) 



where yUe = 2 is the mean mass per electron. Since the torus equi- 
librium temperature stays within the keV range it would take ~ 10' 
years to become non-degenerate, ensuring that the cold torus always 
remains in a degenerate state. 



^ We note that turbulent viscosity is negligible in our case since 
the torus is very dense and metal rich. The corresponding a parameter 
(Shakura & Sunyaev 1973) is very small and is given as cf = v/c^/H = 
6 X 10""'r^|.y/Wio where the disk scale height H is in unist of 10 km. 

We recall that in our model the age of the system is given by fage ~ 
10* yrs for / >> t where t is the magnetic field decay timescale 
given in Niebergal et al. (2006). Here, the period and the spin-down 
rate are units of 10 seconds and 10"'^ s/s, respectively. 



Appendix B: Torus equilibrium temperature 

In general at a given temperature the maximum torus density be- 
low which the matter is non-degenerate, p„^, is found by equat- 
ing the torus temperature to the Fermi temperature, T = Tp^mi = 
123.6 MeVyU^^'Vio^ (see appendix in OLNI), or. 



Pnd - 460;/e g cm ^T^'^'^ 



The corresponding scale height of the torus atmosphere is 



//.,n, = — ~ 12 cm 



A'atm^QS.l^ 



(B.l) 



(B.2) 



where g, = GMqs/( VI/?,)^ ~ 4x 10"MQs.i.4/-Rf ,5 cm s"^ is the eflFec- 
tive gravity at the torus, and 



Vth 



kT 



3.1 X 10' 



cm s-'r"^ 

1/2 -"keV- 



(B.3) 



The corresponding column density of iron atoms in the non-degenrate 
atmosphere is 



» f Ti r 1 n2^ -2 keV t,15 

A'pe = Pnd^/atm ~ 5 X 10^' Cm ^ 

^atm™QS,1.4 



(B.4) 



B.1. Atmosphere opacity 

The torus atmosphere is highly optically thick to radiation (Tatm > 10' 
for temperatures of order 1 keV, where t„ is the total optical depth 
of the atmosphere; e.g. §3.3 in Clayton 1983). For conduction, the 
opacity is lower thus conduction determines the atmospheric temper- 
ature profile. For densities p ~ p„d and temperatures in the keV range 
relevant here, the conductive opacity is found to be of the order of 
/Cc ~ 200 cm^ gm (see §3.4 in Clayton 1983). The resulting con- 
ductive optical depth is = ^CcPndM.aim. ~ 6 X 10^. The standard 
gray atmosphere model gives T^tm t'^* where t is the optical depth 
in the atmosphere. This ensures the effective temperature at the top 
of the atmosphere to be less than 1/10 of that of the bulk body of 
the torus. Since the blackbody cooling rate goes as T'^, the cooling 
is dominated by cooling from the sides of the torus until the area 
of the atmosphere greatly exceed that of the sides of the torus, or, 
Ar, > lO*Rt. Using equation l |A.7t this occurs at a time, in years, 
t > 3.6 X 10^ yrs x R^i^/T^j.^, much later than the time relevant here. 



B.2. Torus thermal evolution 

The torus cooling is given by, 

Lbb, = 47tR;o-T* ~ 2.9 x lO" erg s"' Rl.T^y . (B.5) 

We will assume that the heated torus atmosphere leaks out from 
the sides of height H„d. The accretion rate is then 



R^ y3 

m = 2;r/?,2/f„dP„dVth ~ 1.7 x lO" g s"' ' " 



3/2 

/'aCm^QS.1.4 



resulting on an accretion luminosity of 



Lace = Jlriic ~ 1.5 X 10 erg s 



3/2 

/'atra™QS.1.4 



(B.6) 



(B.7) 



where we estimate that 10% (i.e. 77 = 0.1) of the mass-energy of ac- 
creted matter is released as radiation. 
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The thermal evolution of the torus is given by the heat capacity of 
the torus, 



dT 

dt 



(B.8) 



In the equation above Cy = Mcv is the torus heat capacity with Nt = 
m/Qj^niii) and it's specific heat is c\ = (k^T/ep) x where kg is 

the Boltzmann constant and ep the Fermi energy. The heat capacity of 
the torus is then 



~ 5.2 X 10^' ergK^ 



2/3 



Applying equation l lA.St gives, 
Thus, equation l |B.8t becomes, 

where 



13/6 



'70.1"t.l5 



7.43 



3/2 1/3, J 



(B.9) 



(B.IO) 



(B.U) 



(B.12) 



B 



,2/3 



14.4- 



The equilibrium temperature can be found by setting d/dt = to get 

(B.13) 



T^ec, = - ~0.52keV-5y5 



A'atm^QS.1.4 

The torus reaches the equilibrium temperature almost immediately af- 
ter accretion ensues in timescales of the order of ~ 1/A. 



Appendix C: Wall equilibrium temperature 

Once the wall is penetrated by the magnetic field and detached from 
the torus it will co-rotate with the star. Unlike in the torus itself, where 
the atmosphere is removed from the sides, in the case of the wall the 
atmosphere evaporates directly along the magnetic field lines from the 
ends. The mass-loss rate of the wall through its atmosphere can then 
be written as. 
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m„ = 27r^,Ar„p„_„dV„.,h = 6.85 x 10 g/cc 



t,15 keV 



1/2 r3/2r9/2 
Pw,LUi-n-°s,14 QS 



(C.l) 



Here, 



Pw.nd 



460 g/cc 



-p3/2 

w,keV 



(from Eq. i IBJT i). 



^r„/(;j„ii,n, wjh) is the thermal velocity, and Ar„ = 2nR^m„/p„ due 
to the wall becoming thinner as it loses mass while maintaining a con- 
stant density (as described in § 14.21 ). 

The resulting accretion luminosity L„ = rjm„c^ heats up the wall 
which is meanwhile cooling as a blackbody at a rate of Lw.bb ~ 2.7 1 x 
10^' erg s^^R^i^T^cw assume the wall to cool mostly from its side 
facing the star). The wall's equilibrium temperature is then obtained 
by setting L„ = L„_bb which yields 



r„,eq ~ 4.74 keV 



1/2d3/4 1/2 



1 "t.l5"V-10 



d3/4 n9/4 



/^w,atm"s,14"QS,10 



(C.2) 



